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[ We investigate the complex-temperature singularities of the susceptibility of the 2D Ising 

I model on a square lattice. From an analysis of low-temperature series expansions, we find 
^ ! evidence that as one approaches the point u = Ug = —1 (where u = e~^^) from within the 
complex extensions of the FM or AFM phases, the susceptibility has a divergent singularity of 
the form x ~ + u)~^''' with exponent 7^ = 3/2. The critical amplitude A'^ is calculated. 
Other critical exponents are found to be = = and Ps = 1/4, so that the scaling 
relation a'^ + 2(3s + 7^ = 2 is satisfied. However, using exact results for [3s on the square, 
triangular, and honeycomb lattices, we show that universality is violated at this singularity: 
(3s is lattice-dependent. Finally, from an analysis of spin-spin correlation functions, we 
demonstrate that the correlation length and hence susceptibility are finite as one approaches 
the point u = —1 from within the symmetric phase. This is confirmed by an explicit study 
of high-temperature series expansions. 
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1 Introduction 



Although exact closed-form expressions for the (zero-field) free energy and spontaneous mag- 
netisation of the two-dimensional Ising model were calculated long ago, no such expression 
has ever been found for the susceptibility, and this remains one of the classic unsolved prob- 
lems in statistical mechanics. Any new piece of information on the susceptibility is thus 
of value, especially insofar as it specifies properties which an exact solution must satisfy. 
In particular, it is of interest to understand better the properties of the susceptibihty as 
an analytic function of complex temperature. Several years ago, some results on complex- 
temperature singularities of the susceptibility for the Ising model were reported Here we 
continue the study of complex-temperature singularities of the susceptibility of the 2D Ising 
model. 



2 Generalities and Discussion of Complex Extensions 
of Physical Phases 

We consider the Ising model on a lattice A at a temperature T and external magnetic field 
H defined in standard notation by the partition function 

Z=Y^ e-^H (2.1) 

with the Hamiltonian 

n = -JY.^i^j-HY.^^ (2-2) 

<ij> i 

where cxj = ±1 are the Z2 variables on each site i of the lattice, (3 = {kBT)~^, J is the 
exchange constant, < ij > denote nearest-neighbor pairs, and the magnetic moment /i = 1. 
We shall concentrate here on the square (sq) lattice but also make some comments for the 
triangular (t) and honeycomb (he) lattices. We use the standard notation K = f3J, h = f3H, 

v = tanhK (2.3) 

, = e'^K ^ (2.4) 

l + v ^ ' 

and 

u = z^ = e-^^ (2.5) 

It will also be useful to express certain quantities in terms of the elliptic moduli /c< and 
/c> = l//c<. For the square lattice these are given by 

1 4m 



sinh(2K)2 (1 - uf 



(2.6) 
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and 



We note the symmetries 

K ^ -K =^ {v ^ -V , z^l/z, M— K ^ K} (2.8) 

where = fc< or /c>. The reduced free energy per site is / = PF = — A^~Mimjv^oo In 2' 
(where N denotes the number of sites on the lattice), and the zero-field susceptibility is 
Xo = ^^QH^ |h=0; where M{H) denotes the magnetisation. Henceforth, unless otherwise 
stated, we consider only the case of zero external field and drop the subscript on xo- K is 
convenient to define the related quantity 

X = l3-\ (2.9) 

For the square (sq) lattice, f{K,h = 0) was originally calculated by Onsager 0, and the 
expression for the spontaneous magnetisation M was first reported by Onsager and calculated 
by Yang |Q. Solutions for f{K,h = 0) and M were subsequently given for the triangular 
(t) and honeycomb or hexagonal (he) lattices; for reviews, see Refs. |^, 0. We denote 
the critical coupling separating the symmetric, paramagnetic (PM) high-temperature phase 
from the phase with spontaneously broken Z2 symmetry and ferromagnetic (FM) long-range 
order as K^, and recall that for the square lattice, Vc = z^, = — 1. 

Here we shall study the susceptibility as a function of complex (inverse) temperature, K. 
For our purposes, it is important to discuss generalized notions of phases and thermodynamic 
quantities. We define a complex extension of a phase as an extension, to complex K, of the 
physical phase which exists on a given segment of the real K axis. As noted, e.g., in Ref. 
[|l|], for zero external field, there is an infinite periodicity in complex K under certain shifts 
along the imaginary K axis, as a consequence of the fact that the spin-spin interaction on 
each link < ij > is cXiCXj = ±1. In particular, there is an infinite repetition of phases as 
functions of complex K; this infinitely repeated set of phases is reduced to a single set by 
using the variables v, z and/or u, since these latter variables have very simple properties 
under complex shifts in K: 

K ^ K + niiT =^ {v —>■ V , z —>■ z , u u , kx k^} (2.10) 

ITT 

K ^ K + {2n + l)— ^ {v ^ 1/v , z ^ -z , u^u, k^ k^} (2.11) 

where n is an integer and, as before, kx = /c< or A;>. On a lattice with an even coordination 
number q, it is easily seen that these symmetries imply that the magnetisation and suscep- 
tibility are functions only of u. Because the shift ( p.ll|) leaves u invariant while mapping v 



;i-i;2)2 



(2.7) 
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to it maps a point in the FM phase (and its complex extension) to itself but maps a 
point in the (complex extension of the) PM phase out of this phase. Consequently, when 
studying complex-temperature properties of the model, it is more convenient to start within 
the FM phase, where the various quantities of interest can be expressed as Taylor series in 
the low-temperature expansion variable u. After this study, we shall proceed to investigate 
the properties of the susceptibility in the interior and boundary of the PM phase. It is useful 
to note that a given point vq or zq corresponds, in the complex K plane, to the set 

K = Ko + ni7r (2.12) 

where n & Z and 

Ko = ~\nzo (2.13) 
while a given point uo corresponds to the set 

K = Ko + ^ (2.14) 

reflecting the structure of Riemann sheets of the logarithm. 

The requisite complex extensions of the physical phases can be seen by using the known 
results on the locus of points on which the free energy is non-analytic. For the square lattice, 
these are given by the circles 

v±{e) = ±1 + 2i/2gie ^ ^^(^^^ = ±1 + 2^/2gia; ^2.15) 

for < 6, uj < 2tt. Recall that the property that this locus of points are circles in both the 
high- and low-temperature variables v and z follows because these variables are related by 
the bilinear conformal transformation ( |2.4| ) which maps circles to circles. For later reference, 
these circles are shown in Figs. 1(a) and 1(b), respectively. 

The circles in f or 2; constitute natural boundaries, within which the free energy is 
analytic but across which it cannot be analytically continued. They thus define the complex 
extensions of the physical phases which occur on the real v or z axes in the intervals —Vc < 
V < Vc or Zc < z < 1 / Zc (PM); Vc < v < 1 or < z < Zc (FM); and —l<v< —Vc or 
1/zc < z <oo (AFM). 

Using the general fact that the high-temperature expansions and (for discrete spin models 
such as the Ising model) the low-temperature expansions both have finite radii of convergence, 
we can use standard analytic continuation arguments to establish that not just the free 
energy, but also the magnetisation and susceptibility are analytic functions within each 
of the complex-extended phases. This defines these functions as analytic functions of the 
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respective complex variable {K, v, z, u, or others obtained from these). Of course, these 
functions are, in general, complex away from the physical line —oo < K < oo. 

We shall also need a definition of singularity exponents of a function at a complex singular 
point. In the case of real K, one distinguishes, a priori, the critical exponent which describes 
the singular behaviour at a critical point approached from the symmetric, high-temperature 
phase from the corresponding exponent for the approach from the broken-symmetry, low- 
temperature phase. For a singular point in the complex plane, we shall again distinguish the 
critical exponents describing the singularity as approached from different phases. Thus, for 
the susceptibility x(C)) where ( refers to one of the complex variables listed before, which 
fails to be analytic at one or more singular point(s) {Cs}, if the leading singularity in xiC) 
can be represented in the power-law form 

x(c)s»5. ~ii-c/cr^^- (2.16) 

as ( approaches (s from within the the phase p, we shall refer to (s as a complex singular point 
and 7s,p as the corresponding critical or singularity exponent for the approach to (s from this 
phase. By analogy to standard usage for physical temperature, we shall set 7s, fm = 7s and 
1s,PM = Is to refer to the critical exponents at (s as approached from within the complex 
extensions of the FM and PM phases, respectively. We shall show that for the specific point 
Ms = —1, 7s,AFAf = 7s,FAf- Critical exponents for other quantities at complex-temperature 
singular points are defined in an analogous manner. The locus of points (s where a given 
function is singular in the complex ( plane will not, in general, be a discrete set, in contrast 
to the case for the Ising model on the physical, real K axis. This is illustrated by the locus 
of points ( p.l5|) where the free energy is singular. Even for a function like the magnetisation, 
which, in the complex FM phase where it is nonvanishing, is an algebraic function in ( = 
OT ( = z of the form M = 11^=1 (C ~ Cs)^% the discrete points (g also in general involve 
associated branch cuts, since the exponents /3s are not integers. 

It should be noted that a phase may exist for complex v or z which is not the complex 
extension of any physical phase. An example of this phenomenon occurs in the present case; 
the fourth region, denoted O in Figs. 1(a) and 1(b) constitutes such a phase. 

In contrast to the usual ferromagnetic critical point the Ising model, which can be ap- 
proached only within the PM phase or FM phase (and similarly, the AFM critical point of 
this model, which can be approached only from within the PM or AFM phase), a general 
complex singularity may be approached from within more than two phases. For example, in 
Fig. 1, the singularities at w = ±i, or equivalently, z = ±i, can be approached from within 
the PM, FM, or AFM complex-extended phases, or, indeed, from the region O which is not 
analytically connected to any physical phase. 
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Since for the square lattice x(^) has the symmetry noted above, x(^) = xi~^): it is 
useful to display the complex-extended phases as functions of u. Under the conformal trans- 
formation M = 2^, the circles in Fig. 1(b) are mapped to a single curve, which is a type of 
limagon of Pascal, defined by 

i?e(M) = l + 2^/2cosa; + 2cos2a; (2.17) 

Im{u) = 2^/^ sin 00 + 2 sm2uj (2.18) 

traced out completely for < u < 2ti. In this variable, there are three complex-extended 
phases, as shown in Fig. 1(c): PM, FM, and AFM. The mapping from z to u reduced the 
number of complex-extended phases from the four which are present in the variable z oi v 
to three; since the points in the O phase are related to those in the complex-extended PM 
phase hy z ^ — z, these two phases are mapped to a single phase in the u plane. The points 
z = ±i, at which the PM, FM, AFM, and O phases are all contiguous, are mapped to the 
single point u = Ug = —1 in the u plane. In the terminology of algebraic geometry, the point 
Us = —1 is a singular point, specifically a multiple point of index 2, of the limagon ( p.l8| ) 



forming the natural boundary between the complex phases, whereas all other points on this 
curve, including the PM-FM critical point Uc and the PM-AFM critical point at u = l/u^, 
are regular (ordinary) points of the curve. Here, a multiple point of index n of a curve C is 
a point through which n arcs of C pass (see, e.g., Ref. |^). For later reference, the physical 
critical points are Mc = 3 — 2^/2 = 0.171572875... separating the PM and FM phases and 
= 3 + 2\/2 = 5.8284271... separating the PM and AFM phases. In the complex K plane, 
these correspond to the infinite set of critical points Kc = ±(1/2) ln(l + \/2) +m7r/2, where 
n & Z. Under the transformation u 1/m, the complex-extended PM phase maps onto 
itself, while the complex-extended FM phase maps to the AFM phase, and vice versa. 

Finally, in terms of the elliptic moduli, the natural boundaries have the very simple form 
of the unit circle in the complex or /c> planes: 

fc< = l/ky = e'^ (2.19) 

with < < 27r. These incorporate the symmetries 

u^l/u k.^^ (2.20) 

v^l/v => k^^k^ (2.21) 



where k^ denotes /c<, fc>, or k. Given the inversion symmetry (p.20|) , it follows that the 



transformation ( p.6|) from u to k^ maps both the complex-extended FM and AFM phases 
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onto the same region, which is the interior of the unit circle in the complex /c< plane. The 
complex-extended PM phase is mapped to the exterior of this circle. Under the mapping, 
the actual limagon in the u plane wraps around the unit circle in the /c< plane twice. In 
particular, both the PM-FM critical point Uc^sq and the PM-AFM critical point 1/uc are 
mapped to the single point, A;< = fc> = 1. The complex-temperature singular point Ug = —1 
is mapped to = ky = —1. 

Having discussed these preliminaries, we proceed to study the susceptibility. 

3 Analysis of Low- Temperature Series 
3.1 Analysis of Series for Xr in the Variable u 

The low-temperature series expansion for x for the Ising model on the square lattice is 

oo 

X = 4u\l + Y.CnU^) (3.1.1) 

n=l 

This expansion has a finite radius of convergence and, by analytic continuation from the 
physical low-temperature interval < u < Uc, applies throughout the complex extension of 
the FM phase. Since the factor is known exactly, it is convenient to study the reduced 
(r) function 

oo 

Xr = 2-\-^X = 1 + E CnM" (3.1.2) 
n=l 

The expansion coefficients c„ were calculated to order n = 9 in 1971 by the King's College 
group and were extended to order n = 21 by Baxter and Enting in 1978 |10| (with exact 



coefficients up to n = 19 and nearly exact n = 20 and 21 terms). Very recently, the c„'s have 
been calculated to order n = 2Q (i.e. x to 0(n^'^)) by Briggs, Enting, and Guttmann as 
part of a general calculation of low-temperature series for g-state Potts models with q = 2 
to g = 10 on the square lattice. We have carried out a dlog Pade analysis of this series to 
investigate the singular behaviour of the susceptibility in the complex u plane. (For reviews 
of this method, see Ref. [l^.) As one approaches a complex singular point denoted s on the 
boundary of the complex-extended FM phase from within this phase, x is assumed to have 
the leading singularity (s) 

X{u) ~ A'Jl - u/us\~^'^ (l + ai,,|l - u/us\ + ..) (3.1.3) 



where A'^ and 7^ denote, respectively, the critical amplitude and the corresponding critical 
exponent, and the dots ... represent analytic confluent corrections. One may observe that 
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we have not included non-analytic confluent corrections to the scaling form in eq. (|3.1.3| ). 



The reason is that, although such terms are generally present at critical points in statistical 
mechanical models, previous studies have indicated that they are very weak or absent for 



the usual critical point of the 2D Ising model |jT2|, |13]. The dlog Pade study then directly 
gives Us and 7^. As noted above, the prefactor is known and is analytic, so we actually 
carry out the dlog Pade study on Xr- This study yields evidence for a divergent branch point 
singularity at a particular complex-temperature point, which we denote Us- (The use of Ug 
above referred to generic complex-temperature point (s), of which there might, a priori, be 
more than one; henceforth, we use this symbol to refer to the speciflc point found from the 
Pade study.) The results for Ug and 7^ from the diagonal and near-diagonal approximants 
are listed in Tables 1 and 2, starting with series for Xr to O^u^"^) and going up to 0(m^^). 
We do not flnd evidence for any other complex-temperature singularities within the range 
described by the small-w expansion. From this Pade analysis of the low-temperature series, 
we infer the values 

= -0.998 ±0.002 (3.1.4) 

7; = 1.52 ±0.06 (3.1.5) 

where the uncertainties are estimates. These results suggest the conclusion that as u ap- 
proaches the point 

Us = -1 (3.1.6) 
from within the FM phase, x has a divergent singularity with exponent 

7; = ^ (3.1.7) 

As noted above, this point = — 1 corresponds to the two points Zg = —Vg = ±i; in the 
complex K plane it corresponds to the inflnite set of points given by 

in nin 

^. = -X + — (3.1.8) 

with n E Z. 

Below we shall show that the critical exponent for the inverse correlation length (mass 
gap) describing row or column connected spin-spin correlation functions at this singular 
point is i^'srow = 1- A naive complex-temperature analogue of the usual argument for the 
scaling relation u'{2 — f]) = 7', in conjunction with our inferred value of 7^ in eq. (|3.1.7| ), 
would lead to the further inference that the exponent describing the asymptotic decay of 
the row or column connected 2-spin correlation function at the singular point = — 1 is 
v'srow = 1/2 (where we append the prime to indicate that the calculation of the spin-spin 
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N 


[{N - 2)/N] 


[{N - l)/N] 


[N/N] 


[{N + 1)/N] 


[{N + 2)/N] 


6 




-1.05026 


-1.07608* 


-0.992720 


-0.990763 


7 


-1.09018* 


-1.03635 


-0.990673 


-0.992677* 


-1.00752* 


8 


-1.000(30 


-1.02216 


-1.00363* 


-0.997028 


-0.988513 


9 


-1.00762* 


-1.18678 


-0.981169 


-0.9954185* 


-0.9963225* 


10 


-0.983770 


-1.13876* 


-0.996264* 


-0.995312* 


-0.998649* 


11 


-0.999211 


-1.00398 


-0.998086* 


-0.9977515* 


-0.998400* 


12 


-0.999976* 


-0.997679* 


-0.998044* 


-0.997225* 




13 


-0.999073* 


-0.996210* 









Table 1: Values of from Pade approximants to low-temperature series for Xr starting with 
the series to 0{v}'^). Superscript * indicates that the approximant has one or more nearly 
coincident pole-zero pair(s) closer to the origin than Ug- Our criterion for near coincidence 

is that \Upole — Uzero\ < 10~^. 





[{N - 2)/N] 


[{N-l)/N] 


[N/N] 


[{N + 1)/N] 


[{N + 2)/N] 


6 




1.975 


2.174* 


1.474 


1.455 


7 


2.297* 


1.879 


1.454 


1.474* 


1.631* 


8 


1.563 


1.771 


1.583* 


1.513 


1.419 


9 


1.628* 


1.533 


1.321 


1.497* 


1.508* 


10 


1.358 


1.732* 


1.507* 


1.496* 


1.536* 


11 


1.515 


1.605 


1.528" 


1.521* 


1.533* 


12 


1.554* 


1.523* 


1.528* 


1.517* 




13 


1.543* 


1.503* 









Table 2: Values of 7^ from Pade approximants to low-temperature series for Xr starting with 
the series to 0{u^^). 
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[N/N\ 




1 1 / 




7 


[6/6] 


0.17154017 


1.9 X 10~* 


1.740 


'7 /7' 


U.i / iODUoo 




1 7/1 c; 


[8/8] 


0.17146527* 


6.3 X 10"^ 


1.791* 


[9/9] 


0.17156858* 


2.5 X 10-5 


1.747* 


[10/10] 


0.17157013 


1.6 X 10-5 


1.748 


[11/11] 


0.17157232 


3.2 X 10-6 


1.749 


[12/12] 


0.17157123* 


7.9 X IQ-'' 


1.751* 


exact 


0.171572875.. 





1.750 



Table 3: Values of Uc^ser. a-nd 7' from diagonal Pade approximants to low-temperature series 
for Xr starting with the series to 0{v}^). 

correlation function involves a limit from the complex FM phase). However, we shall show 
that the situation near the complex-temperature singular point « = — 1 is considerably more 
complicated than the case at the physical critical point with its simple scaling relations 
v'{2 — rj) =7' and z/(2 — i]) = 7. Among other things, we shall show that the correlation 
length and x ^t^^ finite when one approaches u = —1 from the complex PM phase. 

The dlog Pade study did not yield any evidence for other complex singularities; i.e., 
it did not give poles whose positions were highly stable as one varied the orders of the 
approximants. As usual, the values of the position of the singular point vary less among the 
Pade entries than the values of the exponent. Also, as expected, the values of 7^ show less 
scatter in the higher-order Pade entries than in the lower order entries. It is true, however, 
that these values of 7^ do exhibit more scatter than the values of the usual susceptibility 
exponent 7' for the PM-FM critical point Uc- To make this comparison quantitative, it is 
sufficient to show the values of Uc and 7' extracted from just the diagonal Pade entries; these 
are given in Table 3. 

3.2 Analysis of Series for g{u) 

Clearly, a property of Pade approximants which is crucial for our study is their sensitivity 
to singularities which are not the closest to the origin of the Taylor series expansion. In 
order to explore the possibility of obtaining a more sensitive probe of the complex singular 
point, we have also carried out a similar study of a series with the physical singularity 
removed. In order to keep the coefficients rational, we actually multiply x by the factor 
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AT 
iV 


I / AT o\/7\rl 


W AT T \ / ATI 

[{N - \)IN\ 


\N/N\ 


\ f AT 1 i\/Arl 

\[N + l)/N\ 


[ / AT 1 o\/7\rl 

[{N + 2)/N\ 


6 




-1.05279 


-1.07822* 


-0.991486 


-1.001045 


7 
( 


1 DQnno* 
— i.uyuuz 


— i.U40UZ 




— u.yoyDo4 


— i.UUooo 


8 


-1.010665 


-1.02141 


-1.00384* 


-0.996939 


-0.989857 


9 


-1.01188* 


-1.34302 


-0.985198 


-0.995437 


-0.996976 


10 


-0.988462 


-1.26478* 


-0.996803 


-0.995048* 


-0.998644* 


11 


-1.00156 


-1.00354 


-0.998084 


-0.997731 


-0.998407* 


12 


-1.001685* 


-0.997653 


-0.998042* 


-0.996954* 




13 


-0.999656* 


-0.994374* 









Table 4: Values of Us from Fade approximants to ^'(m) = (1 — 6m + u'^Y^^Xr, starting with 
the series for g{u) to order 0(m^^). 

[{u - Uc){u - l/u^Y^'^ = (1 - 6m + u'^y^'^ and thus study 

g{u) = {l-Qu + uy/^Xr{u) (3.2.1) 

This is an old technique (see, e.g., Ref. [^). Note that the spurious finite branch-point 
singularity introduced at the PM-AFM critical point has no effect on our analysis, since the 
small- |m I series for x and g only apply in the FM region, which is not contiguous with the 
region of the FM-AFM critical point in the u plane. The results of our Fade analysis of the 
series for g{u) are given in Tables 4 and 5. As one can see from these tables, the resultant 
values for Ug and 7^ are in very good agreement with those from our analysis of the series for 
Xir- One improvement that occurs is a slight reduction of the number of entries with nearly 
coincident pole- zero pairs, which should increase the accuracy of the results somewhat. 

3.3 Analysis of Series for x the Variable A;< 

Finally, it is useful to transform the series for the susceptibility from the usual low-temperature 
variable u to the elliptic modulus fc< and to study the resultant series. An important moti- 
vation for this is that x is given formally by a sum over all connected correlation functions, 
and these correlation functions, which can be computed exactly in terms of certain Toeplitz 
determinants |l^, 0, have explicit forms which are polynomials in the complete elliptic 
integrals K{kx) and E{kx) fl^, where = fc< in the FM and AFM phases and k^ = fc> in 
the FM phase. The variables A;< and /c> are thus natural ones for low- and high-temperature 
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AT 
iV 


[{N - 2)/N\ 


If AT T \ / ATI 

[{N - 1)/N\ 


[N/N\ 


I / AT 1 i\/7\rl 

[{N+l)/N\ 


[ / AT 1 o\/Arl 


6 




1.998 


2.179* 


1.464 


1.561 


7 
( 


o 078* 


i.you 


i .04U 


1 /1/1Q* 


i.04U 


8 


1.664 


1.769 


1.585* 


1.512 


1.434 


9 


1.674* 


1.356 


1.372 


1.497 


1.516 


10 


1.418 


1.490* 


1.513 


1.493* 


1.536* 


11 


1.575 


1.600 


1.528 


1.524 


1.533* 


12 


1.576* 


1.523 


1.528* 


1.514* 




13 


1.551* 


1.481* 









Table 5: Values of 7^ from Fade approximants to g{u) = (1 — 6m + u'^Y^^Xr, starting with 
the series for g{u) to order 0(m^^). 

series expansions of Xy respectively. We therefore have transformed the known small-|M| 
series to one in k^, which takes the form 

00 

X = 2''kl{l + j:c'nk<) (3.3.1) 

n=l 

The series in parentheses defines a reduced function Xr = 4:k~'^x before. Since as a function 
of M, fc< has the expansion near u = Ug = —1 

k< = -l + ^{u + lf + 0{iu + lf) (3.3.2) 

with no linear term, it follows that in the variable A;<, the singular form of x corresponding 
to eq. ( p. 1.31) , as A;< approaches the point fc< = —1 from within the FM or AFM phase (i.e. 



from within the interior of the unit circle in the complex /c< plane) is 

X(A:<) ~5:|l + A;<|~^^/2|^l+6i|l + A;<| + ...) (3.3.3) 

where B'^ is the critical amplitude for this expression of the singularity in terms of the variable 
/c<. We have performed a dlog Fade analysis of the series in fc< for Xr, i-e., an analysis of the 
function d\nXr/dk^. Our results for the diagonal entries are given in Table 6. The results 
from the series in fc< agree very well with those which we obtained from the other series. 
The analysis of the series in /c< also gives results for the regular FM-FM critical point at 
/c< = 1 and the associated exponent 7' of comparable accuracy to that of the series in u. 
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[N/N\ 




ift, \ 111 

1 yl^<)s,ser. ~r i | 


% 


[b/bj 


— u.yys / DO 


1 1 n— 3 

i.Z X iU 


i.Oio 


[7/7] 


-0.999940 


0.60 X 10-4 


1.552 


[8/8] 


-1.000012 


1.2 X 10"^ 


1.555 


[9/9] 


-0.999734 


2.7 X 10-4 


1.543 


[10/10] 


-0.999819* 


1.8 X 10-4 


1.546* 


[11/11] 


-0.999873 


1.3 X 10-4 


1.549 


[12/12] 


c.p. 







Table 6: Values of (A;<)s,ser. and 7^ from diagonal Pade approximants to low-temperature 
series for Xr starting with the series to 0{k^). c.p. denotes a complex pair of poles close to 
— 1. Superscript * indicates that the approximant has one or more nearly coincident pole-zero 
pair(s) closer to the origin than (A;<)s — —1. As before, our criterion for near-coincidence is 

that \{k^)poi(. — {k^)zero\ < 10 4. 

We have also carried out an analysis of the series for Xr in the variable u using differential 
approximants (for further details on this method, see section 7 below). This yields results 
close to those in Table 6. 

One important new piece of information can be obtained from our analysis of the series for 
X in the variable A;< near the singular point /c< = (A;<)s = Ug — —1: this is an answer to the 
question of whether the critical exponent is the same when one approaches this point from 
within the interior of the complex-extended FM phase and from within the complex-extended 
AFM phase. Recall from Fig. 1(c) that in the u plane, these approaches are distinct, since 
the complex FM and AFM phases lie on opposite sides of Ug. Since the complex FM phase 
is the one which contains the origin in the u plane, our analysis of the small- ]ii| series could 
only determine the singular behaviour as u approached Ug from within the complex FM 
phase. However, since the complex FM and AFM phases are mapped onto each other in the 
A;< plane, our analysis of the series in this variable shows that the exponent 7^ is the same 
for the approach to Ug = {k^)s — —1 from within both the FM and AFM phases (as has 
been implicit in our notation). 

Since the low-temperature series for x is not usually given in terms of the variable /c<, 
we note that it has an interesting feature. The first few terms are 

, , 13,2 13,3 139,, 139,, 685,6 2739 , 7 51603, o ^,,„, , 
Xr^l + k^ + -kl + -kl + —k'^ + —k'^ + —k'^ + ^kl + -^k\ + 0{k\) (3.3.4) 
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Near the point Us = (fc<)s = —1, the terms up to order O(fc^) of the series exactly cancel 
among each other in a successive pairwise fashion, so that the first nonzero terms in the 
series for Xr start at order 0{k^) (i.e. for x 0{k^)). To say it differently, these first six 
terms can be expressed as (1 + k^){l + (13/2^)A;| + (139/2^)A;^). One can interpret this as 
being a hint of a structure which persists to all orders in the exact susceptibility; with this 
motivation, one can add and subtract terms in higher orders so as to put Xr in the form 

X = (1 + fc<)/i - 2-''klh (3.3.5) 

Computing the functions f{k^)j, j = 1,2 and defining the convenient variable 

y = 2-^k\ (3.3.6) 

we find 

/i = l + 2Qy + 556?/2 + 10960?/3 + 206412/ + 3775480?/5 + 67668304?/*^ + 1194824896?/^ 
+208565759801/8 + 360778731928?/^ + 61950174438561/^° 
+ 105730294168640?/ii + 17952780821083687/^2 ^ 0{y^^) (3.3.7) 

and 

/2 = 1 + 33?/ + 770?/2 + 15650?/3 + 296006?/^ + 53633357/^ + 945043647/^ + 1633461856?/^ 

+278441539642/8 + 4697355452582/^ + 0{y^^) (3.3.8) 

We have performed a Pade analysis of the functions dlnfj/dy. We find strong evidence 
of a singularity in fi of the form fi ~ |1 — as k"^ increases toward 1 from below. 

Combining this with the (1 — fc<) prefactor, it follows that the singularity in x arising from 
the /i term can be written as 

(1 + A;<)/i ~ |1 - A;<|-^/^|l + k^\~^/^ (3.3.9) 

as approaches 1 from below. For example, the [6/6] Pade approximant gives 16y = = 
0.999844 as the position of the singularity, and 1.745 as the exponent. The results from the 
study of /2 are also consistent with this conclusion for the singularity in x- Recalling that 
|1 + k<:\~^^'^ = const x |1 + u\~^^'^ as A;< — * — 1 or equivalently, as m — —1, one sees that the 
singular form ( |3.3.9| ) agrees very nicely with our determination of the complex-temperature 
singularity by analyses of the small-|M| series for Xr{u), giu), and Xr(^<) given above. 
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r AT 1 ATI 

[iV/iVJ 


Us 


D ( Ai \2/3 


[U/ UJ 


—0 993282 


128915 


[7/7] 


-0.9942585 


0.129496 


[8/8] 


-0.994204 


0.129460 


[9/9] 


-0.994088 


0.1293865 


[10/10] 


-0.995080 


0.130178 


[11/11] 


-0.993497 


0.129294 


[12/12] 


c.p. 


c.p. 



Table 7: Values of {A'^)'^^^ from Pade approximants to small-|M| series for {xrY^'^- c.p. 
indicates a complex pair of poles near to = — 1. 

3.4 Critical Amplitude at Us 

In order to calculate the critical amplitude A'^ in the susceptibility as one approaches u = —1 
from within the FM phase, we compute the series for {XrY^'^'''- Since the exact function 
iXrY^'^'" has a simple pole at Ug, one performs the Pade analysis on the series itself instead 
of its logarithmic derivative. The residue at this pole is — ^^(A^^)^/'''^, where A'^^ denotes 
the critical amplitude for Xr- Using our inferred value 7s = 3/2 to calculate the series and 
the value Ug = —1 to extract A'^g, we finally multiply by the prefactor to obtain A'^ = 
Au'lA'^g = —4:A'g. (Alternatively, one could extract (A^^)^/'^^ from the residue by dividing 
by the measured pole position, Ug^ser. from the given Pade rather than the inferred exact 
position, and could use the prefactor '^u'l to get A'^; the differences between the two 
methods are quite small and vanish asymptotically; these differences are incorporated in the 
final uncertainty which is quoted for the critical amplitude.) Our results from the diagonal 
Pade entries are listed in Table 7. From this analysis, we calculate 

a; = 0.186 ±0.001 (3.4.1) 

where the quoted uncertainty is an estimate. 

This value may be compared with the low-temperature critical amplitude in this model 
at the usual PM-FM critical point, Mc, defined by %(«) ~ A'^\l — uju^""^!^ as u ^ Uc from 
below. A'^ was determined first by analysis of low-temperature series expansions [Q and sub- 
sequently to higher accuracy by analytic methods [|l3l Qto be A'^ = 0.068865538... Using our 



^Ref. [ITl actually gives the critical amplitude A'^j, defined by x(T) - A'^tI^ - Tc/T|~'^/''; we have 
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determination of A'^, it follows that A'^/A'^ = 2.701 ±0.015. We note that this is consistent, to 
within the numerical accuracy, with the analytic relation A'^/A'^ = (— InUc)''''^ = 2.6966995... 

4 Singular Behaviour of Other Quantities at i^^ = — 1 
4.1 Specific Heat 

In this section we extract the singular behaviour of the exactly known thermodynamic quan- 
tities at the complex-temperature singular point Ug = —1. We begin with the specific heat. 
It is convenient to consider K~^C, which (in units with = 1) is given by [0 



2(^K{k) - E{k)) - (1 - + k,'K{k 



(4.1.1) 



where the elliptic modulus is 



_ 2 _ 2 _Au^l'^{l-u) 

and its complementary modulus is = (1 — k^)^/^. We use the standard convention that 
the branch cut for the complete elliptic integrals runs from m = fi;^ = ltom = oo along the 
positive real axis in the complex m-plane. From eq. ([4.1.2| ), it follows that as u approaches 
— 1 with Im{u) positive, k, ioo (and if Im{u) < 0, then k, — > —ioo, taking the usual 
convention for the branch cut for the square root), and thus k' —>■ oo. From inspection of 
( [4.1.1|) , it is clear that as m —1, C diverges, with the leading divergence arising from 
the last term, -{4:/it){k'^ / k^)K{k) {A/ti)k' K{k). Using the identity (see e.g., Ref. fg^l) 
k'K{k) = K{iK/K') and the fact that as A 1, A"(A) (1/2) ln(16/(l - A^)), we can 
express the most singular term as (4/7r) ln(4|K|) as m ^ —1. Next, using the fact that near 
u = —1, 1/k, = — (i/8)(l + m)^ + 0((1 + u)^, we find, finally, that the leading divergence in 
C as M — 1 is 

C r^-K^Jn(-^) (4.1.3) 



Taking the value of Kg on the first Reimann sheet in eq. (|3.1.8| ), i.e.. Kg = —in/ 4, this 
becomes 

C~--lnf^— ) (4.1.4) 

Since the elliptic integrals and also the factor n' only depend on k^, the leading singularity 
is the same whether u approaches Ug = —1 with Im{u) positive, negative, or zero, and 



converted this to A' here. 
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also whether the approach occurs from within the complex FM, AFM, or PM phases. The 
logarithmic divergence in C at Ug is evidently of the same type as the divergence at the 
physical PM-FM and PM-AFM critical points. Note, however, that at these latter points, 
K ^ 1 (/t' — i> 0) so that the k,'K{k,) term (which gives the leading divergence at u = — 1) 
vanishes, and the divergence arises instead from the first term in the square brackets of eq. 
( [4.1.1|) , 2K{k). Another obvious difference is that, while the specific heat is required to be 
positive at physical temperatures, it is, in general, complex at complex-temperature points, 
and the critical amplitude at Ug is real but negative. The critical exponents corresponding 
to this logarithmic divergence in C at Us are 

a, = < = (4.1.5) 



4.2 Magnetisation 

Next, we make use of the exactly known expressions for the spontaneous magnetisation M to 
analyse the behaviour as a function of complex temperature. In particular, we shall extract 
the critical exponent (3s at the complex-temperature singular point u = Us = —1- 

For the square lattice, for real temperature, M vanishes for K < Kc^sq (where for clarity 
we restore here the subscript indicating the lattice type) and, for Kc^sq < -ft' < C)0 is given 
by 0] Af = Msg = (1 - A;|)i/^ or, in terms of u, 

(1-m)V2 ^^■^■'> 

What is normally discussed is the vanishing of M at the usual PM-FM critical point, Mc = 3 — 
2-\/2, with exponent (3 = 1/8. However, as we discussed in section |^, the function describing 
the magnetisation for positive temperature can be analytically continued throughout the 
complex extension of the FM phase, up to the boundaries of this phase, which, for the 
square lattice are specified by the limagon ( |2.18| ). Carrying out this analytic continuation. 



one sees two important results: (1) the only point, other than the physical PM-FM critical 
point, where M vanishes continuously, is at u = = — 1. Defining an associated critical 
exponent as 

\M\ ~ const. X |1 - u/usf" (4.2.2) 
as u approaches Us from within the complex-extended FM phase, we find the value 

Ps,sq = I (4.2.3) 

(2) at all other points (i.e., all points except Uc and Ug) along the boundary of the complex 
extension of the FM phase, M vanishes discontinuously. Result (2) follows because if one 
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starts in the physical PM phase, a similar analytic continuation argument shows that M 
vanishes identically all throughout the complex extension of this phase. Inspection of (|4.2.1|) 
shows that the (real and imaginary parts of the) analytic continuation of M are nonzero as 
one approaches the boundary of the complex-extended FM phase from within that phase at 
points other than u = Uc and u = —1. Therefore, M must vanish discontinuously as one 
crosses this boundary from the complex FM to the complex PM phase, as claimed. 

For the complex-extended AFM phase, we use the well-known symmetry which holds on 
loose-packed lattices: under the transformations K — >■ —K and ai rjiai, where rji = 1 
(— 1) for (Tj on the even (odd) sublattice, the Hamiltonian is invariant, while the uniform 
magnetisation M and the staggered magnetisation Mgt interchange their roles. The above 
transformation takes ^ k^, z ^ 1/z, u ^ 1/u. The expression (|4.2.1| ) is invariant and 
thus describes the staggered magnetisation in the physical AFM phase as well as the uniform 
magnetisation in the physical FM phase. As before, one generalizes this to a definition of 
the staggered magnetisation in the complex-extended AFM phase by analytic continuation 
from the physical region — oo < K < —Kc throughout the complex AFM phase, as indicated 
in Fig. 1. One sees that analogues of the two results which we obtained for M can also be 
derived for Mst'- (1) Mst vanishes continuously at two points on the border of the complex- 
extended AFM phase, namely, u = l/uc = 3 + 2^/2, the usual PM-AFM critical point, 
and u = Us = —1, the complex-temperature singular point where M also vanishes; and 
(2) at all other points on the border of the complex-extended AFM phase, Mst vanishes 
discontinuously as one crosses this border into the PM phase. Note that, as is clear from 
Fig. 1, the only point where the complex FM and AFM phases are contiguous is the single 
point u = Us = —I, or equivalently, the two points z = ±i = —v. 

In passing, we note that for both the uniform and staggered magnetisations, the apparent 
divergence at m = 1 plays no role since this point is outside the two respective regions 
(complex FM and AFM phases) where the expression (|4.2.1| ) for these quantities applies. 



4.3 The Behaviour of the Inverse Correlation Length as u ^ —1 
4.3.1 Approach Prom Within Complex (A)FM Phase 

For the Ising model on the square lattice, in the physical low-temperature phase with real 
K in the interval Kc < K < oo, the asymptotic decay of the row (or equivalently, column) 
connected correlation functions is given by 

< c^(o,o)fX(o,n) >conn.~ n-^e-l"!/^^^^'— (4.3.1) 
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where the inverse correlation length (mass gap) is 

^fIltow = Hiv/zY) = In 













+ z) _ 



(4.3.2) 



We now analytically continue this result into the complex extension of the FM phase and 
inquire where the mass gap vanishes. We find that for points within, and on the border of, 
the complex-extended FM phase, the mass gap vanishes for the following set: 

GM,row = for z = {zc, ±i} (4.3.3) 

i.e., the usual PM-FM critical point Zc = \pl — 1 and at the two points Zg = ±i {us = — 1). 
The additional apparent zero at the PM-AFM critical point z = —1/zc is not relevant for 
the complex FM phase because this point lies outside this phase and thus outside the region 
which can be reached by analytic continuation of the original formula ( [4.3.1| ); however, it 



will be relevant for the inverse correlation length defined within the complex-extended AFM 
phase (see below) . We note the somewhat subtle point that the correct analytic continuation 
of the physical, real-K theory to complex K requires that one use row = ln((f/2;)^) as 
given in eq. ( |4.3.2| ) and not ^pMrow = 2\iii{v/z); although these are identical expressions for 
physical K, the latter form would miss the zero in ^plfrow at u = — 1. 

We now extract the critical exponent (s) for this inverse correlation length (mass gap) at 
z = ±z as these points are approached from within the complex-extended FM phase. These 
exponents (which will turn out to be equal) are defined by 

^FM,row ~ const, x \z ^ i\'^±i,row fgr z ±i (4.3.4) 
from within the FM phase. Expanding ^plj about these points gives 

^Fl,row = + 2(±z -l){zT^) + Oiiz t ^f) (4.3.5) 
from which it follows that 

^'i,row = ^-i,row = 1 (4-3.6) 

This motivates the use of a single exponent to describe the singularity at the single point 
Us = —1 corresponding to z = ±i, as approached from within the complex-extended FM 
phase: 

i^'s,row = 1 (4.3.7) 

By the standard argument noted above which shows that on a loose-packed lattice such 
as the square lattice 

r]ir]j < aiaj > {-K)conn. =< (JiCTj > {K)conn. (4.3.8) 
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it follows that the same inverse correlation length ([4.3.2|) describes the asymptotic decay of 
the connected 2-spin correlation functions along a row or column in the AFM phase, i.e., 
formally, ^pMrow = ^AFMrow^ although the same expression is used in different phases. Now 
^ AFM row vanishes at the physical PM-AFM critical point z = —l/zc, i.e., u = 1/uc, and at 
the complex-temperature singular points z = ±i. It follows that the expansion ( [4.3.5|) also 



controls the critical exponent as one approaches the points u = —1 from the AFM phase 
(i.e. from the left in Fig. 1(b)). 

The value of u'^ is evidently the same as the value v' = 1 for the physical PM-FM 
critical point, as approached from within the FM phase. However, one encounters several 
new features at the complex-temperature singular point, which we now discuss. 

One may also extract a correlation length critical exponent from the asymptotic decay 
of the diagonal {d) correlation function, 

< (^(0,0) (^(n,n) >conn.~ n-^e""/^^"''' (4.3.9) 

where the distance In the physical FM phase, 

rAi4 = -2-^/Mn(fc^) (4.3.10) 



-2-1/2 In 



< 

4m 



(4.3.11) 



As with the row (column) spin-spin correlation functions, we may analytically continue eqs. 
( [4.3.9| )-( [1.3.11| ) to apply throughout the complex extension of the FM phase. Although the 



detailed form of ^fm^ is different from that of ^pMrow^ they both have the same complex 
zeroes within the range of this analytic continuation, i.e. the complex FM phase and its 
boundary. It may be recalled that near the physical PM-FM critical point, 

A;^ = 1 + 2(4 + 3V2){u - u,) + 0{{u - m^)^) (4.3.12) 

so that the exponent which describes the vanishing of the inverse correlation length character- 
izing the diagonal correlation function is that same as that for the row or column correlation 
functions. However, the situation is different at m = — 1: using the expansion near u = —1 
(see also eq. (|3.3.2| )) 



kl = l-]^{l + uf + 0((1 + uf) (4.3.13) 



it follows that near m = — 1, 



i-FM4 = '^'""i^ + ^f + 0{{l + uf) (4.3.14) 
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Hence, the correlation length exponent describing the vanishing of the inverse correlation 
length for diagonal correlation functions, at u = = — 1, as approached from within the 
complex FM phase, is 

K4^ag. = 2 (4.3.15) 

i.e., twice the value of the mass gap exponent extracted from the row/column correlation 
functions. This situation is unprecedented for critical exponents at physical critical points. 

4.3.2 Approach from Complex PM Phase 

From an analysis of the asymptotic decays of both the row/column and diagonal correlation 
functions in the complex extension of the symmetric, PM phase, we find that the correlation 
length does not diverge as one approaches the points v = ±i corresponding to the point 
u = —1 from within this phase. This finding is very important, since it implies that the 
susceptibility is finite at m = —1 when this point is approached from within the complex PM 
phase. 

In the physical PM phase, the row (or column) correlation function has the asymptotic 
decay [|I^ 

< ^^(o,o)^(o,n) >~ |n|-i/2e-l"l/«^^^— (4.3.16) 



where 



CpM,row = Hz/v) = In 



'1 



1 -V 



(4.3.17) 



As before, we may analytically continue this throughout the complex-extended PM phase. 
The mass gap ^PMrow vanishes only at the physical PM-FM critical point Vc = \pl — 1. 
(The apparent zero at —Xjvc = —{V2 + 1) is not relevant because this point lies outside the 
complex-extended PM phase where the above analytic continuation is valid.) In particular, 
as one approaches the points v = ±i from within the complex PM phase, ^pM,row lii(~l)) 
so that < (7(0,0) o"(o,n) >~ {—l)'"'\n\^^^'^ as \n\ oo. 

We find the same result for the diagonal correlation function, which, in this complex PM 
phase, has the asymptotic decay 



< ^(o,o)0^(„,„) >~ Inl'^/^e-'/^^^''" (4.3.18) 



where r = 2^/^ Ira I and 



Ui4 = -2-^/'ln(A:>) (4.3.19) 

:i~^ 



-2-/Mnf,^^ (4.3.20) 
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Now near f = ±i, 

= -1- {vTif + 0{vTiY) (4.3.21) 

Hence, although $,pMd vanishes at the physical critical point Vc, it is finite at the points 
V = ±i, where ^p^d = 2~^/^ln(— 1), so that < a"(o,o) o"(n,n) >~ 1)" at these points, 

just as was true of the row and column correlation functions. One thus encounters precisely 
the type of situation that we discussed before in Ref. [jl|, where Re{C,~^) = 0, but Im{^^^) 
is nonzero. Note that the sum 

E^o(-l)~"^"^^^ (where no is an umimportant lower cutoff) 
is finite. Of course, although the correlation length is finite at f = ±i, as approached from 
within the PM phase, it is singular at these points since it is unequal to the value obtained 
as the points are approached from a different direction in the complex f , or u planes. 
The asymptotic decay of the general 2-spin correlation function < o"(o,o)0"(m,ri) > has been 
calculated (using Toeplitz determinant methods) [|T9] ; carrying out an analytic continuation 
of this result from the physical PM phase into its complex extension, we again find that the 
correlation length is finite at f = ±i. Since a divergence in x on the border of the (complex 
extension of the) PM phase can only arise from a divergence in the sum over 2-spin correlation 
functions contributing to Xi the above results constitute an analytic demonstration that the 
susceptibility is finite at the points v = ±i, as approached from within the complex PM 
phase. This is in sharp contrast to the approach from within the complex FM or AFM 
phases, where we have shown that x is divergent. This type of phenomenon is, again, to 
our knowledge, unprecedented in the study of singularities in thermodynamic functions at 
physical critical points. 

4.3.3 A Theorem on x 

From our results in the previous two subsections, using the same reasoning as in Ref. [|1|, 
we can infer the following theorem: 

Theorem. The susceptibility x has at most finite non-analyticities on the natural bound- 
ary curve (circles in v or z, limagon in u) separating the complex-extended PM, FM, and 
AFM phases, apart from the divergent singularities as one approaches the point Uc = 3 — 2v^ 
from within either the complex PM or FM phase and the point = — 1 from within the 
complex FM or AFM phase. 

5 Scaling Relations and Other Critical Exponents at 

Us = -1 
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5.1 ai + 2/3, + 7,; = 2 



Using our result ( |3.1.71 ) for 7^, together with the exponents a'^ and (3s extracted from the 



known exact expressions for C and M in eqs. ( [4.1 .51) , ( [4.2.3|) , we find that the complex 



analogue of the scaling relation (from the low-temperature side) a' + 2/3 + 7' = 2 is satisfied: 

a', + 2/3, + 7: = 2 (5.1.1) 

where the subscript s indicates that this refers to the point = — 1 and the primes indicate 
that the approach to this point is from the complex-extended broken-symmetry phases, FM 
or AFM. To be precise, this relation is satisfied to within the numerical accuracy of our 
determination of 7^ in eq. ( p.l.5| ) and is satisfied exactly if one uses our inference in eq. 



( p.l.7| ) of the exact value of 7,. However, our results in the previous section, in particular, 
the demonstration that x is finite at f = ±i (u = — 1) as approached from within the 
complex PM phase, and hence that 7, < 0, already shows that 

a, + 2/3, + 7, ^2 (5.1.2) 

i.e. the scaling relation for the approach from within the PM phase, is not valid at Ug. We 
do not know of any extension of the arguments for usual exponent relations to complex 
temperature, so it should not be considered a surprise that such relations do not hold at a 
complex-temperature singularity. 



5.2 Hyperscaling Relations 

Since we have shown above that the inverse correlation length is finite at the points v = ±i 
when approached from within the complex PM phase, the corresponding exponent z/, < 0. 
Hence, the hyperscaling relation dv = 2 — a does not hold at these points, as approached 
from within the complex PM phase. Concerning the hyperscaling relation for the approach 
to Us = —1 from within the complex FM or AFM phases, namely, du' = 2 — a', if one used 
the correlation length exponent i's,row = 1 extracted from the row or column correlation 
functions, then this relation would be satisfied. However, the situation is more complicated, 
since, in particular, u'^ ^^^^g = 2^ v'^^^^^. 

6 Violation of Universality at Ug = —1 

Although the scaling and hyperscaling relations discussed above were found to be satis- 
fied, complex-temperature singularities clearly have different properties from physical, real- 
temperature critical points. Among other things, quantities which are real for physical T 



22 



in general become complex for complex T. Furthermore, various positivity relations, such 
as the property that the specific heat C > is not true even when C is real. One should 
therefore be cautious concerning the question of whether a given property associated with a 
physical critical point will apply at a complex-temperature singular point. Indeed, we shall 
now demonstrate a violation of universality at the complex-temperature singular point Ug- 

We recall the meaning of universality as applied to statistical mechanical models not 
involving frustration or competing interactions: the universality class, as specified by the 
critical exponents, depends on (i) the symmetry group G of the Hamiltonian and the related 
space of the order parameter; (ii) the dimensionality of the lattice, but (iii) not on the details 
of the Hamiltonian, such as additional spin-spin couplings (provided that these are invariant 
under G and do not introduce frustration or competing interactions), and (iv) not on the 
lattice type (again, provided that this does not cause frustration). 

We shall now demonstrate, using exact results, that property (iv) is violated at the point 
Ms = — 1. In order to do this, we use the expressions for the spontaneous magnetisation 
on the triangular and honeycomb lattices. These can be written in the same general form 
( [4.2.1 ) for the square lattice, but with elliptic moduli which are different functions of z: 



M= (l-(fc<,A)2)'/' (6.1 
where instead of the relation (p^ ) for k^^gq, one has 



(1 + 3^2)1/2(1 _ ^2)3/2 



(6.2) 



and 



_ Az^l\l - z + z^f/^ 



These apply to the physical FM phases for each lattice, i.e. where < A;< < 1. The explicit 
forms in terms of the usual low-temperature variable u for the triangular and honeycomb 



lattices are thus [21 



and 22 



_ /1 + m\3/8/1-3m\1/8 , , 

(1 + ^2)3/8(1 _ 4^ ^^2)1/8 
= (l-.)3/4(l + .)l/4 (6.5) 

which apply within the respective FM phases on these lattices and vanish elsewhere. (Recall 
that since the honeycomb lattice has an odd coordination number, q = 3, M and x ^i-re not 
invariant under z —z as they are for lattices of even q.) We first note a similar feature of 
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the spontaneous magnetisation on all three lattices: M vanishes continuously at the same 
generic set of points, namely the respective PM-FM critical points on each lattice, and the 
point M = — 1, or equivalently the two points z = ±i. As is well known, the critical exponent 
j3 = 1/8 is the same at the respective physical PM-FM critical points. However, this is 
not true at m = —1: for the square lattice, the critical exponent was extracted above as 
Ps,sq = 1/4, but for the other lattices, 

3 

Ps,t = Ps,hc = Q Ps,sg (6.6) 



Given the fact that M has the same form (3J.) in terms of the (different) elliptic moduli 



fc<^A for the three lattices, and given that Ug = —1 maps to k^ jy = —1 for each of these 
lattices, it follows that M ~ |1 + fc<,A|^ as fc<^A approaches —1 from with the FM phase for 
each case, with the same value p = 1/8. However, this is not the same as the usual meaning 
of universality, since the fc<^A differ as functions of z for each of the three lattices. 

Since early studies of low-temperature series expansions, it has been known that different 
lattice types have different numbers of complex-temperature singular points (see, e.g., Ref. 
|25| and references therein). However, to our knowledge, the obvious violation of univer- 



sality noted above has not been explicitly discussed in the literature. Indeed, in an early 
study of complex-temperature (and complex-activity) properties of the 3D Ising model [p3| . 



from analyses of low-temperature series expansions on the simple cubic (sc), body-centered 
cubic (bcc), and face-centered cubic (fee) lattices, it was found that the numerical evidence 
was consistent with the equality of critical exponents on these three lattices. The critical 
exponents at the complex-temperature singular point u ^ —0.285 in the Ising model on the 



simple cubic lattice were recently determined to higher precision in Ref. ||2J]. It would be 



useful to calculate longer low-temperature series for C, M, and x on the bcc and fee lattices 



to compare with the higher-accuracy critical exponents obtained in Ref. pj] for the sc lat- 
tice. However, our exact results on jSg already show that universality does not, in general, 
hold at complex-temperature critical points. 

Some possible insight into this violation may be gained by remembering that even at 
physical critical points, universality does not, in general, hold when there is frustration. One 
of the earliest examples is the (isotropic) antiferromagnetic Ising model on the triangular 
lattice, for which there is no PM-AFM phase transition at finite K. Accordingly, when 
examining a given singular point to see if one could expect universality to hold, one of 
the first things which one would necessarily check would be the presence or absence of 
frustration, which, in turn, would involve checking whether various spin configurations only 
partially minimize the internal energy. But this initial check cannot be performed in the 
usual way at a complex-temperature singular point, since at such a point the Hamiltonian 
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and internal energy are not, in general, real numbers. 



7 Analysis of High- Temperature Series Expansion for 
the Susceptibihty 

In section 4, as a consequence of our study of the complex-temperature behaviour of corre- 
lation lengths, we showed analytically that the susceptibility is finite at the points v = ±i 
(i.e., u = —1) as approached from within the complex extension of the PM phase. In this 
section we shall carry out a study of high-temperature series expansions for the suscepti- 
bility. The results confirm our analytic demonstration and give further information about 
X at these points. To our knowledge, this is the first time that a comparison has been 
made of the behaviour at a complex-temperature singularity as approached from both the 
complex-extended FM (AFM) and PM phases. For technical reasons, the study of the high- 
temperature series in the vicinity of v = ±i turns out to be considerably more difficult than 
was the case with the low-temperature series in the vicinity of the equivalent single point 
u = —1. We begin with a simple dlog Pade study, which is adequate to confirm the absence 
of a divergent singularity; we then proceed to a study with differential approximants. 
Recall that the high-temperature series expansion for the susceptibihty is given by 

oo 

X = 1 + E (7-1) 
n=l 



in terms of the usual high-temperature expansion variable. We have also transformed the 
series to one in the elliptic modulus variable {k 



vl/2 



>J 5 



X = 1 + E <(^>)"^' (7-2) 

n=l 



via the relation 



Note the symmetry 



{k^y/' = 2v/il-v') (7.3) 



The motivation for using the variable (fc>)^^^ is the same as was discussed in reference 
to the low-temperature series, viz., that the exact expressions for the spin-spin correlation 
functions which actually contribute to the susceptibility are polynomials in the complete 
elliptic integrals of modulus fc> in the PM phase (multiplied by algebraic functions of (A;>)^/^) 
||T6| , |l7|, |T^. Under the mapping from v to (A;>)^/^, the boundaries between the phases 
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transform as follows: the circle v = —1 + is mapped to the right-hand unit semicircle 

in the (fc>)^''^ plane, i.e., (A;>)^/^ = e*"^ with -ir /2 < (f) < it/2. Given the symmetry ([7^) , 
this is a two-fold covering; in particular, the image of the the PM-FM critical point Vc and 
the point —1/vc is the single point (fc>)^/^ = 1. Similarly, the circle v = l-|-2^/^e*^ is mapped 
by a two-fold covering to the left-hand unit semicircle, (k^Y^'^ = e*"^ with (3/2)7r > > 7r/2. 
The points —Vc (the PM-AFM critical point) and 1/vc, are taken to the point (A;>)^^^ = — 1. 
Finally, the points v = ±i which lie on the intersections of the two circles are mapped to 
(k^Y^"^ = ±i, respectively. 

For the square lattice, the a„ have been calculated to the very high order f by Nickel 
26| . We have performed a dlog Fade analysis on this series and have found evidence against 



a divergence in x as f approaches ±i from within the PM phase. Since x is real for real 
V, it follows that if the Fade approximant for dln{x)/dv has a pole at f = f o with residue 
Rq at some complex point Vq, then it also has a pole at v = Vq with residue Rq. Writing 
the singular part of x as (1 — v/vo)~"'° near v = vo and recalling that Rq = —70, it follows 
that at the two complex-conjugate poles, the real parts of the exponents are equal, while the 
imaginary parts (if nonzero) are reversed in sign. Thus, without loss of generality, it suffices 
to consider only the singularity at v = i. 

We find that the Fade approximants to d\n{x)/dv yield a reasonably stable pole near 
to V = i, with Re{'-fv=i) < 0. However, even for rather high order [N/N] Fade entries with 
15 < < 23, the pole position is not as close to the singular point as one would require for 
accurate results; typically, \vi^ser. — ^| — 0.08, much larger than the usual level of O(10~^) (or 
better) which one would expect for reasonable accuracy. 

We have therefore studied the equivalent series ( [7.2[ ) in the elliptic modulus variable, 
(A;>)^/^. In terms of the latter variable, the leading form of the singularities at (A;>)^/^ = ±i 
is given by 

since the Taylor series expansion of (A;>)^/^ =1= z as a function of v, near the points v = ±i, 
starts with the quadratic term: 

(fc>)i/2 = ±z± '-{v T if + 0{{v T if) (7.6) 

Of course, the Fade approximants exhibit the well-known pole at (A;>)^^^ = 1 due to the usual 
FM-FM critical singularity and the sequence of poles and zeroes starting near to (A;>)^/^ = — 1 
and continuing outward along the negative real axis attributed to the finite (1 + x) In |1 + x| 
FM-AFM singularity, where x = (A;>)^/^. For reference, the high order dlog Fade entries for 
the pole position near the physical FM-FM critical point get this accurate to order 0(10^^). 
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I AT / AT^ 

[N/N\ 




\/ u ^ 1/2 • 1 

\{i^>)i -A 


Re{'yi) 


[7/7] 


0.972846z + 0.032631 


4.2 X 10~^ 


-0.4297 


[8/8] 


1.00313/ + 0.031169 


3.5 X 10"- 


-0.5749 


[9/9] 


1.006322 + 0.015344 


1.7 X 10"^ 


-0.7720 


[10/10] 


1.00729i + 0.003183 


0.80 X 10"^ 


-0.83235 


[11/11] 


0.997106i + 0.001684 


3.3 X 10"-'^ 


-0.6024 


[12/12] 


0.999979i + 0.000388 


3.9 X 10"^ 


-0.6432 


[13/13] 


0.999715i + 0.000650* 


7.1 X 10~^ 


-0.6404* 


[14/14] 


1.00053z + 0.006735 


6.8 X 10~^ 


-0.7166 


[15/15] 


1.00026i + 0.006752* 


6.8 X 10"-'' 


-0.7101* 


[16/16] 


0.998362i + 0.003564 


3.9 X 10"-'^ 


-0.6380 


[17/17] 


0.998506Z + 0.003520* 


3.8 X 10"^ 


-0.6409* 


[18/ i8J 


u.yyyo/uz + u.uuouuy 


Q n v 1 n— 3 

o.U X iU 


— U.D / Z / 


[19/19] 


0.998932i + 0.002885 


3.1 X 10-3 


-0.6442 


[20/20] 


0.999033i + 0.002739* 


2.9 X 10-=^ 


-0.6449* 


[21/21] 


0.999021i + 0.002416 


2.6 X 10-3 


-0.6391 


[22/22] 


0.999472i + 0.002638 


2.7 X 10-3 


-0.6567 


[23/23] 


0.999487i + 0.002646 


2.7 X 10-3 


-0.6573 


[24/24] 


0.999162i + 0.002051 


2.2 X 10-3 


-0.6374 



Table 8: Values of (A;>)/ and Re{ji) from diagonal dlog Padc approximants to high- 
temperature series for x starting with the series to 0{{k]l'^Y^) = 0{v^^). 

They also yield extremely precise determinations of the exponent, 7; indeed, among the 
[N/N] entries with N around 20, the values of 7 only differ from 1.75 by amounts of order 
10^^ — 10^^. As regards complex-temperature singularities, the approximants exhibit two 
poles which converge to (A;>)^^^ — ±i These, together with the values of Re{'ji), are shown 
in Table 8. The values of Re{'^i) are stable and are negative, indicating that x has finite 
singularities at = ±i as these points are approached from within the complex PM phase. 
(The values of /m(7i) will be discussed below.) This confirms the conclusion of our analytic 
study in section 4, viz., that x is finite at v = ±i as approached from within the complex 
PM phase. (That it is singular is obvious since it has different values when approached from 
different directions in the complex z, or u plane.) 

To check the sensitivity of the dlog Pade analysis to possible divergent singularities, we 
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have studied the test function 



X{x)t,st = ^,.(1 - x)-'/^ - AafmA'^ + x) ln(l + x) + + x^Y^^'^ + B (7.7) 

where x = (A;>)^/^. This function incorporates the known PM-FM and PM-AFM singulari- 
ties, a hypothetical divergent singularity at (A;>)^/^ = v = ±i with exponent corresponding 
to (1 + v'^)~'^% and an additional background term B. We have used the known critical 
amplitudes Ac,x = 2^/*^^^^, = 1.4153665..., where A 

c,v is defined by \sing ~ 
for V Vc and is given by Ae,, = 0.77173406... [0, |7|; and Aafm^x = 2^/^A,^afm = 0.28, 
where A^^afm is defined by Xsing ~ — ^i),afm(1 + ^^/'^c) ln|l + v/vc\ for — 1> y —Vc and is 
given by A^^afm — 0.20 We have varied A^ a. and B over a range of values and 7^ over 
the range 1/4 < 7<j < 3/2, and have found that if these quantities have values such that 
the dlog Fade approximants locate the singularities at f = (/c>)^''^ = ±^ with an accuracy 
comparable to that which we observe in Table 8, then the approximants also yield reason- 
ably accurate values for 7^. In particular, if we make the values of Ag^x, Is and/or B so 
small that the Fade fails to yield an accurate value for 7^, then it also fails to locate the 
singularities at (/c>)^^^ = ±i in the test function with the accuracy that it does successfully 
locate them for the actual x- Hence, the dlog Fade would not miss divergent singularities 
at f = (kyY^'^ = ±i if they were really present in 

However, the dlog Fade method is not, in general, satisfactory for a finite singularity, 
and hence, given that -Re(7i) < 0, one cannot, a priori, trust the actual values of 7^ which it 
yields. The appropriate technique to investigate such finite singularities in the presence of 
background terms is provided by differential approximants (DA) |jl5|, ^ ^ |3l|, |32| . To 
be precise, if a function of a generic variable ( is of the product form /(C) = YijLii^—C/CjY^^ ? 
then the dlog Fade method is, in principle, adequate to obtain the positions and exponents 
for the singularities even if some pj < 0. If / has an additive background term near a 
singularity, i.e., if fsing ~ ^(1 — CKj)~^^ + -^(C) with B{Q analytic near then, if a given 
Pj > 0, the first Darboux theorem [|14| shows that for sufficiently high orders, this will 
dominate over B{() so that the dlog Fade method can still be satisfactory, but if the given 
Pj < 0, then one should use differential approximants. 

We recall that in the differential approximant method the function f = x being approx- 
imated satisfies a linear ordinary differential equation (ODE) of K'th order, £M,L/if(C) = 
j:f=oQjiC)D'fK{C) = RiO, where Qj{0 = E^oQi/^' and i?(C) = EloReC 0, S 



p2| . In one implementation of the method |3^, |3T]], D = d/dQ, while in another |32, 15 
D = Cd/d(. We adopt the choice for D used in Ref. [0, |T^; these authors have found that 
both choices give comparable results. The solution to this ODE, with the initial condition 
/(O) = 1, is the resultant approximant, labelled as [L/Mq; ...-jMk]- The general solution of 
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the ODE has the form fj{C) ~ Aj {QIC - 01"^' + ^(C) for C 0- The singular points Cj 
are determined as the zeroes of Qk{0 and are regular singular points of the ODE, and the 
exponents are given by —pj = K — 1 — <5/<-i(Cj)/(CiQx(Ci))- Further details on the method 



can be found in Refs. [|T^, |3T|, For an extrapolation procedure to be discussed below, 
we shall use a number of poles at different positions close to the singularity; for this reason, 
we use unbiased differential approximants. Studying the susceptibility series in the variable 
V, we find that the differential approximants do not yield singularities sufficiently close to 
±i to be accurate, just as was true of the dlog Fade method (which is a special case of DA). 
As before, we have obtained considerably better results with the series in the elliptic mod- 
ulus variable (fc>)^/^. We have calculated the K = 1 differential approximants [L/Mo]Mi] 
for 4 < L < 24 and 10 < Mq < 20 with Mi = Mq, Mq ± 1 subject to the constraint 
L + Mo + Ml + 1 < 49 (terms up to 0{v^^) = 0{{k]/'^y^) were used). Many of the poles 
may reflect finite singularities along the arcs of the circles bounding the complex FM phase, 
as discussed previously [|I|]. To consider a pole to represent the singularity at (k^Y^'^ = i, 
we require that its distance from this point satisfy |(A;>)^/^ — i\ < 1 x 10^^. Secondly, we 
require that the pole position lie within the circle |fc>| = 1 which forms the boundary of 
the complex FM phase, since the high-temperature series only converges to x within this 
boundary. In Table 9 we display the K = 1 differential approximants with L even which 
satisfy these conditions. (The approximants with odd L are not listed to save space; they 
yield conclusions in agreement with those obtained from the approximants with even L.) 
Evidently, there is a large scatter of values of 7^. Also, the distances from the singularity 
are usually larger than the O(10~^) level which one would normally consider necessary for 
accurate results. However, we can still draw useful information from this table. First, all of 
the values of 7^ satisfying the two conditions above have the property that Re^'ji) < 0. Sec- 
ond, the approximants in Table 9 which yield poles closest to (A;>)^/^ = i, namely [10/16; 18], 
[20/12; 12], [22/10; 12], and [24/10; 11], do give roughly consistent values of Re{-fi). If we 
plot the values of Re^'ji) as a function of the distance |(A;>)^/^ — i\ and extrapolate to zero 
distance from the singularity, we obtain Re{'^i) ^ —0.65. This is consistent with the values 
obtained for this quantity from the dlog Fade study. 

We can make further progress by noticing an important correlation: the sign of /m(7j) is 

1 /2 

opposite to the sign of the deviation i?e((A;>)j ) from zero. Indeed, when we plot the values 
of /m(7j) as a function of Re{{ky)i), they can be roughly fit to a line going through zero when 
Re{{ky)y'^) = 0. But the singularity which we are studying is at (/c>)^/^ = i, so we are led 
to the tentative conclusion that for this singularity, Im{'-).j) = 0. It follows also that for the 
conjugate singularity at (/c>)^/^ = —i, Im{'j_i) = 0, whence 7^ = 7„j = 7^. As regards the 
values of Im{'yi) from the dlog Fade study, one can see that all of the high-order diagonal 
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[L/Mo;Mr] 




\{k>)T-i\ 


li 


[4/18; 


19J 


0.993373^ + 0.004293 


0.79 X 10^2 


-1.420 - 1.079i 


[4/18; 


20J 


0.991349i + 0.002545 


0.90 X 10^2 


-1.820 - 0.8796i 


To / 1 o 

[8/18; 


16J 


0.990874i + 0.003398 


0.97 X 10"^ 


-1.950 - 0.3293i 


[8/18; 


17] 


0.988550i - 0.006664 


1.3 X 10"^ 


-1.944 + 1.741i 


[8/18; 


18] 


0.990839i + 0.004293 


1.0 X 10~^ 


-2.164 - 0.3726i 


[10/16 


;18] 


0.999122i - 0.001421 


1.7 X 10"^ 


-0.8568 + 0.4512i 


[12/10 


;10] 


0.9987905i - 0.005881 


0.60 X 10"^ 


-1.243 + 0.4854i 


[14/10 


;io] 


0.990038i + 0.007285 


1.2 X 10"^ 


-2.294 - 1.181i 


[14/12 


;14] 


0.990566i - 0.001845 


0.96 X 10"^ 


-1.611 + 0.3934i 


[18/10 


;io] 


0.990270Z + 0.002944 


1.0 X 10-2 


-1.907 - 0.3510i 


[18/12 


;12] 


0.998597i + 0.01127 


1.1 X 10-2 


-0.7024- 1.417i 


[20/10 




0.994749i - 0.001108 


0.54 X 10~2 


-1.267 + 0.3876i 


[20/10 


;12] 


0.995176i - 0.002265 


0.53 X 10-2 


-1.214 + 0.5585i 


[20/12 


;10] 


0.998042i + 0.008232 


0.85 X 10-2 


-0.1837- 1.433i 


[20/12 




0.995593i - 0.003949 


0.59 X 10-2 


-1.105 + 0.9183i 


[20/12 


;i2] 


0.998991i - 0.000353 


1.1 X 10-^ 


-0.6833 + 0.3278i 


[22/10 


:11] 


0.997787/ - 0.002213 


3.1 X lO-'' 


-0.81 19 + 0.6317/ 


[22/10 


;i2] 


0.9982122- 0.001353 


2.2 X lO-''^ 


-0.7283 + 0.5223^ 


[24/10;9] 


0.993917i- 0.006415 


0.88 X 10-2 


-1.850 + 0.9510i 


[24/10;10] 


0.997266i - 0.001653 


3.2 X 10-3 


-0.9898 + 0.4162i 


[24/10;ll] 


0.997902i + 0.0004882 


2.2 X 10-3 


-0.5841 - 0.0597i 


[24/10;12] 


0.997133i - 0.008122 


0.86 X 10-2 


-1.120 + 1.803i 



Table 9: Values of {k>)i and 7j from K = \ differential approximants to high-temperature 
series for %. See text for definition of [L/Mo,Mi] approximant. 
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approximants are characterized by the same (positive) sign for Re{{k-^Y/'^ ^ i.e. the pole 
positions are shghtly to the right of (fc^)^/^ = i. This correlates with the observed feature 
that these Pade approximants have the stable nonzero negative values Jm(7j) ^ —0.25. 
Indeed, when we examine other poles in the dlog Pade approximants reasonably close to the 
point {ky.y/'^ = i, we observe the same correlation that Im{'ji) has a sign opposite to that 
of i?e((A;>)^/^). This suggests that if we had a reasonably large set of such nearby poles, 
then we could carry out an extrapolation similar to the one that we performed with the 
differential approximants. There are not enough close poles with i?e((A;>)^/^) < to do this 
accurately, but one can say that such an extrapolation is crudely consistent with /m(7j) = 0. 

Summarizing, then, the results of our analysis with differential approximants, like those 
of the simple dlog Pade study, confirm our analytic demonstration that x is finite (although, 
of course, singular) at the points v = (kyY^"^ = ±i (i.e., u = —1) when approached from 
within the complex PM phase. Furthermore, if we restrict to the four poles closest to 
V = {k^Y^"^ = i, the DA's yield values for Re{'yi) which are roughly mutually consistent. 
Note that for these four poles, the distance from the singularity, 2 — 3 x 10"'^, is comparable to 
that for the poles produced by the high-order dlog Pade approximants in Table 8. Moreover, 
when we extrapolate the values of -Re(7j) to apply precisely at the location of the singularity, 
we obtain a value consistent with that from the dlog Pade approximants. Using a similar 
extrapolation procedure with the differential approximants, we infer that Jm(7j) = 0. 

We recall that in principle, x have singularities elsewhere on the arcs forming the 
natural boundary of the complex PM phase, as discussed in Ref. JIl. It was proved there that 
if such singularities exist, they must be finite. Indeed, from a simple dlog Pade study, it was 
observed in that paper that the poles lay close to these arcs. The differential approximants 
presumably give some information about these possible finite singularities along the arcs. 
We shall leave the more detailed investigation of the behaviour on these arcs to future work. 
However, from our present results, we may make some interesting observations. To the 
extent that the poles in the differential approximants away from (fc>)^'^^ = ±i lie near to 
the arcs bounding the complex PM region, they might reflect finite singularities in x along 
these arcs. The results of our study then suggest that the associated exponents would have 
nonzero imaginary parts, at points along the arcs in v or circle in (fc>)^^^ apart from the 
points (A;>)^^^ = ±1 (v = ±Vc) and (k^Y^"^ = ±i (v = In assessing this possibility, 

it is, of course, incumbent upon one first to check whether there are any rigorous theorems 
forbidding this. We have not been able to derive any such theorem. Of course, the usual 
rigorous theorems governing the behaviour of thermodynamic quantities and, in particular, 
their critical exponents, assume physical values of the temperature. Two properties which 
must be satisfied are that x must be real and positive for physical temperature. These two 
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properties do not exclude the existence of a complex exponent at a complex-temperature 
singularity in x- To demonstrate this, we consider a generic form for the singular part of x 
with such singularities: 

Xs^nAC) = ^-(1 - C/Co)-^" + Alii - C/Co)""° (7.8) 

where ( denotes a generic variable {v, (/c>)^/^ z, etc.), Imi^o) ^ 0, and the form (fTSp applies, 
say, for |^| < |^o|- Evidently, this has the property that 

Xsmg(C ) Xsmg(C) (7-9) 

which implies that if is real, then x((^) is also real.0 We have also checked that one can 
write down examples of Xsing of the form in eq. (|7]B|) which yield not just real, but positive, 
X for physical temperature. If x has a singularity at the point C = Co of the leading form in 
eq. ( |7.8| ) where the real and imaginary parts of 70 are given by 

70 = 7^+^7/ (7-10) 

then for C = Co(l + ee^^), as e ^ 0, 



X 



cos(7, ln(|e|)) - % sin(7, ln(|e|)) 



(7.11) 



Given our theorem above that that x must be finite on the arcs away from v = and f = ±i 
(approached from within the FM or AFM phases), it follows that 7^ must be negative away 
from these points. The effect of the nonzero 7^ is to produce an infinitely rapid oscillation 
in X as C ^ Co with |C| < |Co|- Although the points on the arcs are not isolated singularities, 
but instead, points on the natural boundaries between the complex-extended PM, FM, and 
AFM phases, this behaviour is quite analogous to the infinitely rapid oscillations at isolated 
essential singularities such as in the functions Csiii(l/C) or Ccos(l/C) at C = 0. Given our 
results that x diverges at u = — 1 as one approaches this point from the complex FM or 
AFM phases, but has a finite singularity as one approaches the point from the complex 
PM phase, one may observe that this behaviour is somewhat reminiscent of the function 
exp(l/C), which has an (isolated) essential singularity at C = and diverges (vanishes) when 
this point is approached from the positive (negative) real axis. 

A physical argument for the divergence in the susceptibility at a second-order phase 
transition, such as the PM-FM critical point, is to note that the magnetisation vanishes 



^ Actually, real C may still correspond to complex temperature; for exa mple, as eq. ( 2.1C ) shows, K = 
r + niTT with r real corresponds to real u and hence real x- The symmetries ( 2.10| ) and ( 2.11 ) imply that for 



the indicated infinite sets of complex values of K, quantities such as x are still real. 
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continuously as one approaches this point from within the FM phase, so that, in the hmit, 
an arbitrarily small external field H has an arbitrarily large relative effect on the resultant 
magnetisation M[H). This argument also agrees with the divergence in the susceptibility 
that we found at the complex-temperature point u — —1, as approached from within the 
complex-extended FM phase, since this is the only other point at which M vanishes con- 
tinuously, starting from within the complex FM phase. For the physical PM-FM critical 
point the same argument motivates the divergence in x one approaches this point from 
within the FM phase. However, our analytic demonstration (and confirmation via high- 
temperature scries analysis) that x has a finite singularity when one approaches u = —1 
from within the complex PM phase shows that for complex-temperature singularities, this 
reasoning for physical temperatures is not applicable in a naive manner. 

8 Conclusions 

In summary, we have carried out a study of the complex-temperature singularities of the 
susceptibility of the 2D Ising model on a square lattice. From an analysis of low-temperature 
series expansions, we have found evidence that as one approaches the point u — Ug — —1 from 
within the complex extensions of the physical ferromagnetic or antiferromagnetic phases, x 
has a divergent singularity. Our results are consistent with the conclusion that the critical 
exponent for this singularity is 7^ = 3/2. We have also calculated the critical amplitude. 
However, from an analysis of the asymptotic decays of spin-spin correlation functions, we 
have shown that the correlation length is finite, and hence the susceptibility is finite (although 
both are singular) at the points v = ±i corresponding to w = — 1 when approached from 
within the complex-extended symmetric, paramagnetic phase. This is confirmed by a study 
of high-temperature series expansions. 

Our results are of interest because they elucidate the behaviour of the susceptibility as 
an analytic function. The goal of calculating the susceptibility of the 2D Ising model (or 
even making a conjecture for this function which agrees with available series expansions) 
has remained elusive for half a century since Onsagcr's calculation of the free energy. Our 
results should be useful for this quest because they provide new properties which must be 
satisfied by such a conjecture or calculation. 

One can think of a number of further related studies to perform. In one direction, using 
low- and high-temperature expansions, we have carried out analyses with dlog Fade and 
differential approximants to investigate complex-temperature singularities of x iii the Ising 
model on the triangular and honeycomb lattices. Our results will be reported in a sequel 
paper. One could also examine other models like the d — 2, q — ?, Potts model. Yet 
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another topic is how ideas of conformal field theory, which have greatly illuminated the 
critical behaviour at the physical critical points of various 2D models, can be applied to this 
complex-temperature singularity. Clearly, there is much work for the future. 
This research was supported in part by the NSF grant PHY-93-09888. 
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Figure Caption 



Fig. 1. Phases and associated boundaries in the complex variables (a) f , (b) and 
(c) u, as defined in eqs. (|27^ )- (|2.5| ). In the variable /c< defined in (|2.6|) , the complex FM 
and AFM phases are mapped into the interior of the unit circle, and the FM phase to the 
exterior of this circle. 
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